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ABSTRACT 


The stationary transport load acts on the surface of the cavity or on the inner surface of the shell reinforcing cavity. The speed of the 
load movement is assumed to be subsonic, which corresponds to the modern transport speeds in the underground facilities under 
study. To describe the motion of a half-space and thick-walled shells, dynamic equations of the theory of elasticity in Lamé 
potentials are used, and for thin-walled shells, the classical equations of the theory of thin shells are used. Equations are recorded in 
the moving coordinate system associated with the load. 
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1. INTRODUCTION 


In the theoretical aspect, the solution was based on the works [1,2,3]. In [4], the first and second boundary-value problems of the 


theory of elasticity for a half-plane with a point source of stationary waves concentrated within it, the potential of which is 
represented in terms of cylindrical functions, are solved by the method of expanding potentials on plane waves. And in [5], using this 
approach, the problem of the stationary load on the contour of a circular hole in a half-space was solved. Using the idea of these 
papers on the superposition of solutions and the re-expansion of plane waves into series in cylindrical functions, in [6], in contrast to 
the exact analytical solution for the subsonic case, when the velocity of the moving load is less than the velocity of the Rayleigh 
waves. 


Statement of the problem for a circular tunnel 

Using the model approach for research, we will represent the tunnel as an infinitely long circular cylindrical cavity with a radius r = R, 
located in a linear viscoelastic, homogeneous and isotropic half-space x < h (Figure 1) parallel to its horizontal boundary (the earth's 
surface). We define the reaction of a half-space on a moving with a constant subsonic velocity c along the cavity surface in the 
direction of the load axis Z P. 


Figure 1 The calculated scheme of a reinforced tunnel 


For this, we use the equations of motion of an elastic medium in vector form [9, 10] 
2- 
~W2- ma ~ is = O Uu 
UV U+(A + D)graddiFu = Pie , (1) 
t 


Here uu, ’ U,, ’ U. ) - vector of displacement of points of the medium; p - material density; u — displacement components; 
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j- Poisson's ratio; 
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E- Operational modulus of elasticity, which have the form [11,12]: 


ple) = F,| 0)—[Re(t—r)pl0de | a 
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y(t)- arbitrary time function; R; (t- t)- relaxation core; Eo - instantaneous modulus of elasticity; We assume the integral 


—1OR 


terms in (5) to be small, then the functions g(t) = (te 7 where y(t)- a slowly varying function of time, @, - real constant. 


Further, applying the freezing procedure [13], we note relations (2) as approximations of the form 


Eg=E\l-L°(o,)-iT* (a, Jp, 


where |" “(x)= J R(e)cos o,tdrt, VT * (@g)= J R(z)sin Opt dT, respectively, the cosine and sine Fourier images of the 
0 0 


relaxation core of the material. As an example of a viscoelastic material, we take three parametric relaxation nuclei 
= —Bt pyl-a ; : rr oe: 
R(t) = Ae /t . On the influence function R(t = T) the usual requirements of inerrability, continuity (except for € = 


T), sign of uncertainty and monotony: 


dR(t r 
R>0O, mo <0, o< f R@dr<1. 
0 


u - the vector of displacements of the environment. 


Since the steady-state process is considered, the strain pattern is stationary with respect to the moving load. Therefore, it is 
convenient to move to a moving coordinate system n = z— ct, connected with the load P. 
Then equation (1) can be rewritten in the form 


1 1 1 ‘ 
Mw. M2 grad divas pVu= (3) 
P s Ss 


2. TASKS OF THE ACTION OF MOBILE LOADS ON AN UNREINFORCED TUNNEL 


In the theoretical aspect, the solution was based on the papers [15, 16]. In [14], the first and second boundary-value problems of the 
theory of elasticity for a half-plane with a point source of stationary waves concentrated within it, the potential of which is 
represented in terms of cylindrical functions, are solved by the method of expanding potentials on plane waves. And in [15], using 
this approach, the problem of the stationary load on the contour of a circular hole in a half-space was solved. Using the idea of 
these papers on the superposition of solutions and the re-expansion of plane waves into series in cylindrical functions, in [16], in 
contrast to the exact analytical solution for the subsonic case, when the velocity of a moving load is less than the velocity of the 
Rayleigh waves. 


3. STATEMENT OF THE PROBLEM FOR A CIRCULAR TUNNEL 


Using the model approach for research, we will represent the tunnel as an infinitely long circular cylindrical cavity with a radius r = R, 
located in a linear viscoelastic, homogeneous and isotropic half-space x < h (Figure 1) parallel to its horizontal boundary (the earth's 
surface). We define the reaction of a half-space on a moving with a constant subsonic velocity c along the surface of the cavity in the 
direction of the Z-axis of the load P. 

Since the steady-state process is considered, the strain pattern is stationary with respect to the moving load. Therefore, it is 
convenient to move to the mobile coordinate system n = z-—ct, connected with the load P. 
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Then equation (1) can be rewritten in the form 


Here Mp = c/Cp, Ms = c/cs - Mach numbers; oe wa +2u Vp ,C,= VL/p — complex propagation velocities of expansion 


waves - compression and shear in a medium. 


When the load acts on the cavity surface, we have 


r= P(8,n), jJ=r,9,n, (5) 


om 


where o,;— components of the stress tensor in a medium, P;(8,n) — components of the intensity of the mobile load P(8,n). 
Since the boundary of the half-space is free from loads, x = h 


O,, =9,, =0,, =0. (6) 

We transform equation (1) by expressing the displacement vector of an elastic medium through Lame potentials 
u = grad ~, + roty. (7) 
Potential y can be represented in the form [7] 
w= Qe, + rot(p;e,, ), (8) 


where ey — ort axis 1. 


With this in mind, (5) takes the form 
u = grad div @, + rot(9,e, )+ rot rot(e,e,, ) (9) 


It follows from (3) and (8) that the potentials @; satisfy the modified wave equations 


00, 
Vo, =M; ie! (10) 


J J On? 


Here M, = Mp, M2 = M3 = Ms. 
We express the components of the stress-strain state (VAT) of the medium through the potentials qj. 


The components of the vector u (7) in cylindrical (8) and Cartesian (9) coordinate systems [17]: 
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a9, 109, A 
u, = ®; i M> 3 . 
Or r 00 Onor 
2 
pa 0p, 1 OP; (11) 
r 00 Or rd0onod 
‘i _ OY, me 0°, ; 
| On s On? " 
— oO”, 4, 4: 0°Q, 
* @x Oy = axdn’ 
pes OY, OP, 4 0’o, 
*  @y Ox OyOn” 
é @* 
u, = an? pas 
on on 
where m2 =|- M? 
Volumetric strain 
e=divu=V’Q,. (12) 


Using Hooke's law, taking into account (9), (11), we find expressions for the stress tensor components in cylindrical and 
Cartesian coordinates 
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Thus, to determine the components of the stress-strain state of the medium, it is necessary to solve equations (9) together 
with the boundary conditions. 


In the moving coordinate system, we apply to the equations of motion and the boundary conditions a complex Fourier transform of 
the form [16] 


PE)=[ Pla) "dn (14) 


oln)=—5— J ole e" ae 


Writing general solutions of the transformed equations of motion of the tunnel in the form (4) - (14), we find the following system of 
algebraic equations for determining the dimensionless trans formants of displacements of an intermediate surface 


— EU, +1€G,W, =-¢* Gio; (15) 


= 2 4 
iEG|U , ~A.G367w, + [1+ “e wo — 


a 1—G, CG, 
2k y 


Wy = C2Fo; 
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where G, =G,/G,;k =h/a; P, =P,al Eh, 


3 
= -Co=Cc) 2} 
Wo Wo} = Ue Whi 0 (34 


The stress at the boundary of the soft layer and elastic among (r = b) in the dimensionless form has the form: 


_ Orm co sntl —_ a 
o,. -{= eee!) sin + yb @ ®) on nOre dé 
a ron An 
. 2 i Ba’ 
t= f2- : 


Ba? H (Ba) + anf Fe Ho (fa) — BaH,,( pe) 


-S|a+me, (@ a)-a aH, (a@ a)cos6|- 


wo gat - nH” (@ a)t+(@ a)H,_,,(a@ a)| 
os SE See le 


7 2), An 


n=2 


osn@ Seis” 

2de 6, =—4nH\ (@ a)H, (B a)+(1+n)a@ aH, (@ a)+H, (PB a), 
An=n@ aH, (@ a)H,,(B a)+n(B a)H,_(B a)H, (@ a)- 
~apaH" (a@ aH” (B a) 


n-l n-l 


Here 5 = p/ps is the ratio of the density of the environment to the density of the soft layer; @, Bp - are functions of € and n. 


We find the following expression for the load transformer, which is transferred to the shell from the side of the soft layer 


4A B ; 
4j|k4= 
Cy = i] K_1=0°73 j Lc yee pen Bay 
y / a= 
“i det||A,, det //A,, | 
Elements of the determinant det|A,, is computed then formula A,, =—2M,,; Ay =—-a)15 A; = nM»; 


Aa = =Ais: Az, = =§,Aji5 A,» = A), *ky(z, \k,(z,); 
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1 
Ags = Ais * Ig (Zp )/Ta(z2)s Aog = Ais *kq(z, 0k, (2) Ay =5A Ay = an 


Ag, = 1 ko(z3)/k,(z4)-2A,, (z,/M,); A, =A,/n; A, =-A,/0;; 
Ay = mo(z,)/T,(z,)—2M,S,(z,)/1,(z, } 
Ajyy = -2M},(ky(z;)/ ky (ze) + 1, (z5 )/ L,(Z, /(z, 1M): 
Aggy = —2M 3 (Iq (zs )/ ky (25) +1 (26d, (26 (ze 1M): 
me m= Jf1-M?; mm, =i-M2; 
Zz, = M,N; Z = M1; z, =Myn; 


<4 =my(l+k,,) f5 =), ky = (b—a)/a; 


Ky, k, - modified Neumann functions; I, [, - modified Bessel functions; The general solution of the equations of the 


motion of the environment has the form (c, at Oe Cs) 
g(r, E) = A, (é)k,, (m,ér)+ DB. (EN, (m,ér) 


y(r, é) =G, (Ek, (m,,ér) + D, (EI, (m,,ér) (16) 


The expression for the trans formant of normal displacement has the form 


(17) 


i=o 


wets | f° Ailacoslen)- cnt 


m “* a’ +f° |det|A,,| 


Define A, (7=12 deseo 5) is obtained from det|A,, 


by replacing j = 20 by the column C with the elements {0; 0; 1; 0; O}. 


After this function A(é) disaiee D(é) from (23) can be calculated from formulas 


2 1 1 A! 1 
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A} = 2M yy + Pi)malG i(k = 1,.......4) 


J 


m,, — minors of the element  Aje. For a specific value of the load velocity C, the denominators under the integral expressions in 


ie 
formulas (14) are transcendental functions with respect to ¢ C real coefficients depending on C, as well as on the mechanical 
parameters of the shell and the layer. Analysis of the integrals of treatment must begin with consideration of cases 

D(E,C, ) =0, which is equivalent to the construction of the dispersion relation in the corresponding problem of propagation of 
free waves and the determination of the denominator from the dispersion curves of the roots for the chosen velocity of the load C. 
at C< Cs are possible for cases. Figure 2 shows the change in the movement of the filler, depending on the thickness of the bodies 


for different values of the rigidity of the aggregate. As can be seen from the drawing (y = 100, 50,10,2), that for a sufficiently rigid 
layer (y = 100), the deflections of the shell essentially decrease. 


Figure 2 Shell deflections as a function of thickness 


For a given speed C, there are one or two different denominator roots. 

2. For some values of C, the denominator has a double root. This case corresponds to a minimum of the corresponding dispersion 
curve in Fig. Such a velocity is called resonance and is denoted by C*% A resonance effect appears, or which deflections and 
contact pressures tend to infinity. 

3. For a given value of C, the denominator has no roots on the real axis, as seen in Figure 2, this will be either, C<Cy (up to 
resonance mode). At this speed of motion, the inversion integrals are not special and can be found by effective numerical 
methods. 
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Dividing the integral (21) into two terms 


1 po 1 pe 
Wy) =— [, x (Q)dQ o Wo= —| x, (Q)dQ 18) 
fa Yaa | 


The value of the integral (18) was found numerically, using the Rmberg method [2]. When the integral is calculated by the Romberg 


method, it is necessary to repeatedly calculate the integrand function. The inverse Fourier transform (29) was numerically fulfilled. It 


is shown that at an integration step of 1.01, the error of the procedure does not exceed 0.3-0.5%. 


4. CONCLUSION 


1. To describe the behavior of viscoelastic materials with unstable properties that do not obey the principle of temperature-time 


analogy, the non-singular singular kernel of heredity is proposed. 


2. A universal algorithm for solving the problems is proposed. 
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